Optimal Error Correction in Topological Subsystem Codes 
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A promising approach to overcome decoherence in quantum computing schemes is to perform active quantum 
error correction using topology. Topological subsystem codes incorporate both the benefits of topological and 
subsystem codes, allowing for error syndrome recovery with only 2-local measurements in a two-dimensional 
array of qubits. We study the error threshold for topological subsystem color codes under very general external 
noise conditions. By transforming the problem into a classical disordered spin model, we estimate using Monte 
Carlo simulations that topological subsystem codes have an optimal error tolerance of 5.5(2)%. This means 
there is ample space for improvement in existing error-correcting algorithms that typically find a threshold of 
approximately 2%. 
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Quantum computing promises to fundamentally further the 
bounds of computability, particularly in such fields as com- 
plexity theory and cryptography, and, in particular, the simula- 
tion of chemical and physical systems. Unfortunately, imple- 
mentations of quantum computing proposals require precise 
manipulations of quantum systems which are highly suscepti- 
ble to external noise. The technical feasibility of any quantum 
computer design thus heavily relies on efficient quantum error 
detection and recovery. This can be achieved, for example, 
by redundantly encoding quantum information in a code sub- 
space of many physical qubits |[TJ— [3] . Such a suitable subspace 
is defined in terms of stabilizer operators JH [H — products of 
individual Pauli operators — and their corresponding eigenval- 
ues. 

Because stabilizers need to be measured during the error re- 
covery procedure, geometric locality of the involved qubits is 
essential for practicality. Topological error correcting codes 
Hi achieve this by arranging qubits on a topologically 
nontrivial manifold with stabilizers acting only on neighbor- 
ing qubits. These codes promise a reliable approach to quan- 
tum computing, because of their stability to errors [I12l - ll8ll : 
A sizable fraction of physical qubits needs to fail before the 
logical information encoded in the system is lost beyond error 
correction. 

To determine the error stability of topologically protected 
quantum computing proposals it is customary to map the 
error correction procedure onto the thermodynamic behav- 
ior of a disordered classical (statistical-mechanical) spin sys- 
tem 

(mum. 

There is a fruitful synergy between quan- 
tum computation and statistical mechanics: On the one hand, 
the stability of quantum computing proposals can be studied 
with the well-established machinery from statistical physics 
of complex systems, and on the other hand, it also opens the 
door to exotic applications of statistical models. 

Unfortunately, there is one caveat: The stabilizers for sur- 
face codes (such as the Kitaev code |@]) and topological color 
codes [7] involve multiple qubits — four in the case of the Ki- 




FIG. 1 : (Color online) Graphical representation of the qubit arrange- 
ment for topological subsystem color codes on a regular triangular 
lattice. Each of the triangular unit cells (large gray triangles) con- 
tains three physical qubits (red balls). The two-qubit gauge genera- 
tors a w a w are shown in green (w = x), yellow (w = y) and blue 
(w = z). These are the lines connecting the qubits (red balls). They 
are arranged such that each physical qubit has two generators of z 
type, one of x type and one of y type. See main text for details. 



taev code, six or eight for color codes. This immensely com- 
plicates physical realizations. However, in stabilizer subsys- 
tem codes |2(1 [2]]] some of the encoded logical qubits are 
"gauge qubits" where no information is encoded. This pro- 
vides ancilla qubits to absorb decoherence effects and, in par- 
ticular, allows breaking up the required measurements for er- 
ror recovery into several individual measurement that involve 
a smaller number of qubits lEoiElll e.g., two. Hence, physi- 
cal realizations are more feasible at the price of requiring ad- 
ditional qubits. Note that extensions and variants have also 
been proposed 

A true advantage is given by topological subsystem codes 
which combine the robustness of topologically based im- 
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FIG. 2: (Color online) (a) A regular triangular lattice satisfies the 
vertex three-colorability requirement (indicated by A, B, C). (b) To 
construct a topological subsystem code, we place three qubits (red 
balls) inside each of the triangular unit cells and connect them with 
<j z <J Z gauge generators (dotted blue lines). The links between 
these triangles are assigned a x cr x and a y a y gauge generators 
(yellow and green solid lines, respectively), (c) For the mapping, 
gauge generators represented by colored lines in (b) are associated 
with Ising spins s x,y,z and the qubits with interactions, (d) Introduc- 
ing new Ising spin variables s zz — s z s' z allows for the removal of 
local Z2 symmetries. 



plementations with the simplicity of subsystem codes where 
only measurements of neighboring qubits are required for re- 
covery. As in the case of surface and color codes, the ideal 
error stability for topological subsystem codes can be com- 
puted by mapping the error recovery problem onto a classi- 
cal statistical-mechanical Ising spin system where the disorder 
corresponds to faulty physical qubits. Here, using large-scale 
Monte Carlo simulations we compute the ideal error correc- 
tion threshold for topological subsystem color codes affected 
by depolarizing noise. Our results show error correction is 
feasible up to 5.5(2)% faulty physical qubits. Remarkably, 
existing error correcting algorithms only reach a threshold of 
approximately 2% I24II25IL leaving ample room for improve- 
ment. 

Topological subsystem codes and mapping. — A stabilizer 
subsystem code is defined by its gauge group Q. Its ele- 
ments are Pauli operators that, by definition, do not affect en- 
coded states. Namely, two states p and p' are equivalent if 
p = Y2i9iP9i w i m 9i an d 9i elements in the algebra gener- 
ated by Q. 

Topological subsystem color codes (9|] are constructed by 
starting from a two-dimensional lattice with triangular faces 
and three-colorable vertices. Here we consider the triangular 
lattice shown in Figs. [T] and Oa). As indicated in Figs. Q] and 
[2b), there are three physical qubits per triangle and the gauge 
group has 2-local generators Gi of the form a w <g> a w , where 
w = x, y, and z. 

Any family of topological codes shows a finite threshold for 



a given local noise source. In other words, when the intensity 
of the noise is below the threshold, we can correct errors with 
any desired accuracy at the price of choosing a large enough 
code in the family. We are interested in the error threshold of 
topological subsystem codes under the effects of depolarizing 
noise, where each qubit is affected by a channel of the form 



V p (p) = (l-p)p- 



E 

w=x,y, 



cr^pa" 



(1) 



Here p represents the density matrix describing the quantum 
state of the qubit and p G [0, 1] its the probability for an error 
to occur. The depolarizing channel plays a fundamental role 
in quantum information protocols where the effects of noise 
need to be considered, e.g., in quantum cryptography [|26il27 1. 
quantum distillation of entanglement [28], and quantum tele- 
portation Q. 

It is expected that there exists a threshold value p = p c such 
that in the limit of large codes, for p < p c error correction suc- 
ceeds with probability 1 and for p > p c the result is entirely 
random. Remarkably, for topological codes in general, one 
can relate p c to a phase transition in a suitably-chosen classi- 
cal disordered Ising spin model, as we detail next. 

To construct the related classical statistical-mechanical sys- 
tem, we place an Ising spin = ±1 for each gauge generator 
Gi . Single qubit Pauli operators a w are mapped onto interac- 
tion terms according to the generators Gi with which they do 
not commute, giving rise to a Hamiltonian of the general form 



H T (S) 



- J E E rru 



9ij 



(2) 



j w—x,y,z 



Here i enumerates all Ising spins and j all physical qubit sites, 
respectively. For each spin Si the exponent G {0, 1} is 
[1] if [anti] commutes with Gi. The signs of the couplings 
t™ = ±1 are then quenched random variables satisfying the 
constraint t/tJt^ = 1. For each j, they are all positive with 
probability 1—p and the other three configurations have prob- 
ability p/3 each. 

In our specific case the Hamiltonian has the geometry de- 
picted in Fig. [2c) and thus takes the form 



(3) 



where j enumerates qubit sites and spins are labeled, for each 
j, as shown in Fig. [2 Notice that z-labeled spins are arranged 
in triangles, and that flipping each of these triads of spins to- 
gether does not change the energy of the system. Therefore, 
there is a Z 2 gauge symmetry. We fix the Z 2 gauge symmetry 
and at the same time simplify the Hamiltonian by introduc- 
ing new Ising variables Sj z = SjSj. Notice that these spins 
are constrained: If j, k, I are three-qubit sites in a triangle, 
s*fsl z s 

E 



ZZaZ k z s zz = 1. The simulated Hamiltonian therefore reads 



T y .S V .S ZZ 



■T Z S X S V 



(4) 
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Note that the Hamiltonian in Eq. has no local symmetries, 
but a global Z 2 x Z 2 symmetry. Indeed, we can color spins 
according to their nearest colored vertex in the original lattice 
[Fig.Oa)], producing three sublattices A, B, and C. Flipping 
the spins of two of these sublattices together leaves the energy 
invariant, giving rise to the indicated symmetry. 

We are thus left with a random spin system with two param- 
eters, T and p. It is expected that for low T and p the system 
will be magnetically ordered. In the ground states each sub- 
lattice has aligned spins and thus the sublattice magnetization 
is a good order parameter: 



— V 



(5) 



where N-p = L 2 /3 (L the linear system size) represents the 
number of spins in one of the sublattices. The threshold p c 
for topological subsystem codes is recovered as the critical p 
along the Nishimori line Bill 



4/3 J = In 



p/3 



(6) 



where the ferromagnetic phase of a sublattice is lost 

Numerical details. — We investigate the critical behavior 
of the classical Ising spin model [Eq. ©] via large-scale par- 
allel tempering Monte Carlo simulations Jill |33D. Both spin 
states and interaction terms are bit encoded to allow for effi- 
cient local updates via bit masking. Detecting the transition 
temperature T c (p) for different fixed amounts of disorder al- 
lows us to pinpoint the phase boundary in the p - T phase di- 
agram (Fig. [4]). 

We choose periodic boundary conditions keeping in mind 
the colorability requirements. Then we can use the magnetiza- 
tion defined in Eq. (0) to construct the wave-vector-dependent 
magnetic susceptibility 



Xm(k) 




ik Hi 




(7) 



where (■ • -) T denotes a thermal average and is the spatial 
location of the spin Si. From Eq. © we construct the two- 
point finite- size correlation function, 



2sin(fc min /2) V [Xm(k 



(8) 



where [• • ] av denotes an average over disorder and k m j n = 
(2-k/L, 0) is the smallest non-zero wave vector. Near the tran- 
sition £l is expected to scale as 



£ L /L ~ X\I}I V {T - T c )\ 



(9) 



where X is a dimensionless scaling function. Because at the 
transition temperature T = T c , the argument of Eq. © is zero 
(up to scaling corrections) and hence independent of L, we 
expect lines of different system sizes to cross at this point. If, 
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FIG. 3: (Color online) Crossing of the correlation function £l/L 
with a disorder rate of p = 0.048. The data exhibit a clear crossing at 
a transition temperature of T c (p) ~ 1.251(8) H3QTI . The shaded area 
corresponds to the error bar in the estimate of T c (p) . Note that error 
bars are calculated using a bootstrap analysis of 500 resamplings. 
Corrections to scaling are minimal at this disorder rate, but increase 
closer to the error threshold. 



TABLE I: Simulation parameters: p is the error rate for the depo- 
larizing channel, L is the linear system size, A^ sa is the number of 
disorder samples, t eq = 2 b is the number of equilibration sweeps, 
T m in [T max ] is the lowest [highest] temperature, and iVr the number 
of temperatures used. 
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however, the lines do not meet, we know that no transition 
occurs in the studied temperature range. 

When determining the transition temperature T c (p) for a 
given disorder rate p, the correlation functions £l/L are ob- 
tained by averaging over several disorder realizations (gov- 
erned by p) for every system size L. Because we are only able 
to investigate limited system sizes L < oo, a careful analysis 
of finite- size effects is required when estimating the transition 
temperature in the thermodynamic limit. 

In all simulations, equilibration is tested using a base-2 log- 
arithmic binning of the data: Once the data for all observables 
agree for three logarithmically sized bins within error bars we 
deem the Monte Carlo simulation for that system size to be in 
thermal equilibrium. The simulation parameters can be found 
in TableU 

Results. — For the pure system (p = 0) there is a sharp 
transition visible directly in the sublattice magnetization. The 
transition temperature T CjPure w 1.65(1) has not been com- 
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FIG. 4: (Color online) Computed phase diagram for the classical 
disordered spin model shown in Eq. (|3}. Each data point T c (p) on the 
phase boundary (dashed curve separating white and shaded regions) 
is calculated by locating the crossing in correlation function £l/L 
for different system sizes L at a fixed disorder rate p. The Nishimori 
line (blue solid line) indicates where the requirement for the mapping 
[Eq. (|6])] holds. The error threshold p c ~ 0.055(2) is found where 
the Nishimori line intersects the phase boundary between the ordered 
phase (shaded) and the disordered phase (not shaded, larger T andp). 
Below p c « 0.055(2) error correction is feasible. The (red) shaded 
vertical bar corresponds to the statistical error estimate for p c . 

puted before. For larger amounts of disorder, a possible transi- 
tion can be located precisely by means of the two-point finite- 
size correlation function [Eq. (0)]. Sample data for a disorder 
strength of p = 0.048 (i.e., this would mean that on average 
4.8% of the physical qubits have failed) are shown in Fig. [3] 
indicating a transition temperature of T c (p) = 1.251(8). At 
p = 0.055(2), the lines only touch marginally such that both 
the scenario of a crossing as well as no transition are compat- 
ible within error bars. For error rates p > p c , the lines do not 
meet, indicating that there is no transition in the temperature 
range studied. 

The crossing of the critical phase boundary T c (p) with the 
Nishimori line [Eq. ©] determines the error threshold to de- 
polarization. Our (conservative) estimate is p c « 0.055(2). 
Our results are summarized in Fig. |4] which shows the esti- 
mated phase diagram. 

Summary. — We have calculated numerically the error re- 
silience of topological subsystem codes to the depolarizing 
channel by mapping the error correction procedure onto a 
statistical-mechanical Ising spin model with disorder. The 
large critical error rate of p c = 5.5(2)%, combined with a 
streamlined error recovery procedure that requires only two- 
qubit interactions, constitutes a promising implementation 
concept for quantum computing. 
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